Complex sine-Gordon Equation 
in Coherent Optical Pulse Propagation 



< 



Q-Han Park* and H. J. Shin* 

Department of Physics 
and 

Research Institute of Basic Sciences 
Kyunghee University, Seoul, 130-701, Korea 

It is shown that the McCall-Hahn theory of self-induced transparency in coherent optical pulse 
propagation can be identified with the complex sine-Gordon theory in the sharp line limit. We 
reformulate the theory in terms of the deformed gauged Wess-Zumino-Witten sigma model and 
■ address various new aspects of self-induced transparency. 
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' Self-induced transparency (SIT), a phenomenon of anomalously low energy loss in coherent optical pulse propagation, 
t— I | was first discovered by McCall and Hahn 0] and the integrability of the SIT equation was demonstrated by employing 
5_! ■ the inverse scattering method Q . When phase variation is ignored in the case for a symmetric frequency distribution 
Qh| ff(Aw) of inhomogeneous broadening, McCall and Hahn have proved an area theorem for pulse propagation. In the 
sharp line limit where the pulse spectral width is much broader than the atomic line width such that g(Aw) = 
5(w — wo), the SIT equation reduces to the well-known sine-Gordon equation and the 2tt area pulse becomes a 1- 
soliton of the sine-Gordon theory. However, when phase variation is included, the area theorem no longer holds and 
the structure of SIT in general has not been well understood except for the construction of explicit solutions by the 
^ , inverse scattering method Q ||. In view of various physical situations where phase variations become important, e.g. 
— ■ chirping effect or the slow reshaping of pulses in the off resonant case it is important to undestand the analytic 
properties of SIT with phase variation. In addition, despite its integrability, the SIT theory in the sharp line limit has 
not been identified with a known 1+1 dimensional integrable field theory, which made a systematic understanding of 
SIT in terms of a lagrangian field theory impossible. 

The purpose of this Letter is to show that the SIT theory with phase variation can be identified with the complex 
sine-Gordon theory in the sharp line limit. The complex sine-Gordon theory, a generalization of the sine-Gordon 
theory with a phase degree of freedom, can be reformulated in terms of a nonlinear sigma model which is known 
as the integrably deformed gauged Wess-Zumino- Witten(WZW) model associated with the coset SU(2)/U(1) [| [§. 
This allows us to address various new aspects of SIT in terms of characteristics of the complex sine-Gordon theory; e.g. 
topological v.s. non-topological solitons, local gauge symmetry, the U(l)-charge conservation, the chiral symmetry 
and the Krammers-Wannier duality for dark v.s. bright solitons. In physical terms, the local gauge transformation 
accounts for the off-resonance effect of SIT as well as the inhomogeneous broadening of medium. The U(l)-charge 
characterizes the effect of phase variation in optical pulses. We also address how our effective potetial formulation 
can be extended to more physical cases without the sharp line limit. 
The SIT equation is given by 

ft : 

dE + 2/3<P> =0 
3D - E*P - EP* = 
dP + 2iAwP + 2ED = (1) 

where Aw = w — wo , d = d/dz , 8 = d/dz ,z = t — x/c,z = x/c. E,P and D represent the electric field, the 
polarization and the population inversion respectively. The bracket denotes an averaging over the distribution function 
of inhomogeneous broadening. Since Eq.(l) is invariant under the interchange (/3,E,P,D) <-> (—/3,E,—P,—D), we 
assume the coupling constant (3 to be positive which we set to one by rescaling E, P and D. In a simpler case 
where phase variation is ignored to make E real and the frequency distribution is sharply peaked at the carrier 
frequency(Aw = 0), we may parametrize E, P and D by 

E = E*=dip , <P>=P = -sin 2<p , L> = cos2<^. (2) 
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Then the SIT equation reduces to the well-known sine-Gordon equation 



dd(p-2/3shx2(p = 0. (3) 

In order to include phase variation as well as the off-resonance effect (Aw =/= 0), we assume E to be complex and 
require that the distribution is sharply peaked not necessarily at the carrier frequency, i.e. g(Aw) — S(Aw — £) for 
some constant £. Introduce a more general parameterization of E, P and D in terms of three scalar fields ip, 9 and rj, 

E = e i(e-2 V )( 2d ™!W_ id } p _ ie Ke-2 n ) sin2 D = C os2^. (4) 
sm if 

The main result is that the SIT equation given in Eq.(l) then changes into a couple of second order nonlinear 
differential equations known as the complex sine-Gordon equation, 

ddip + A^^dridr] - 2 sin 2tp = (5) 
sin (p 

ddn : — — (drjdip + dridip) = (6) 

sm 2ip 

together with a couple of first order constraint equations, 

2 cos 2 ipd-q - sin 2 ipdO - 2£ sin 2 ip = (7) 
2 cos 2 tpdri + sin 2 ipdO = . (8) 

The complex sine-Gordon equation first appeared in 1976 in a description of relativistic vortices in a superfluid 
, and also independently in a treatment of 0(4) nonlinear sigma model m. Note that Eqs.(5) and (6) consistently 
reduce to the sine-Gordon equation when phase variation is ignored so that rj = 0, 6 = ? an d the system is on 
resonance(£ = 0). Earlier works on the complex sine-Gordon theory have focused only on Eqs.(5) and (6). However, 
the constraints in Eqs.(7) and (8) which are new expressions constitute an essential part of the SIT theory, particularly 
in connection with a local U(l)-gauge symmetry of SIT as explained later. Thus we will call Eqs.(5)-(8) as the complex 
sine-Gordon theory. The gauge symmetry structure as well as the integrability of Eqs.(5)-(8) may be best understood 
if we reformulate the complex sine-Gordon theory in terms of the action principle. The proper action is given by a 
group theoretical nonlinear sigma model action known as the deformed gauged WZW action defined as follows; 

S(g, A, A, 13) = S wzw {g) + ^(-Adgg- 1 + Ag^dg + AgAg' 1 - AA) - 5 potential 

^potential = ^ / ^gTg- lr T (9) 

where Swzw(g) is the conventional SU(2)-WZW action and g is an SU(2) matrix function. Tr denotes the trace and 
T = —T = i<73 = diag(i, — i) for Pauli matrices Oi. The local gauge fields A — a(z, z)as, A = a(z, z)oz are introduced 
to gauge the U(l) subgroup of SU(2). Owing to the absence of the kinetic terms, the gauge fields A, A act as Lagrange 
multipliers which result in the constraint equations. One of the nice properties of our fomulation is that the equation 
of motion arising from the action (9) takes a zero curvature form, 

S g S = -i- J Tr[ d + g- x dg + g~ x Ag + (3 XT , 8 + A + ^Tg }g- 1 Sg = (10) 

where the constant A is a spectral parameter and the square bracket denotes the commutation. The constraint 
equations coming from the A, A- variations are 

S A S=±- f Tr( -Bgg- 1 + gAg' 1 - A )5A = (11) 
2w J 

5 A S=^- { Tr( g- x dg + g~ x Ag - A )SA = . (12) 

27T J 

The action (9) is known to possess the local U(l)-vector gauge symmetry under the transform; 

g^bT x gh , A^ A + h~ x dh , A^ A + h~ x dh (13) 
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where h = exp(^(z, 2)03), as well as the global U(l)-axial vector gauge symmetry under g — > /15ft, for a constant ft. 
In order to identify Eqs.(10)-(12) with the SIT equation, we fix the vector gauge by choosing 

A = £T , A = (14) 

for a constant £. Such a gauge fixing is possible due to the flatness of A, A ||. Also, we parameterize the 2x2 matrix 
9 by 

/ 2ir] ■ • id \ 

= e iri<7 3 giip(cos6<7 1 —sin8a 2 ) e iri<T3 _ I e COS tSlTlipe \ (15) 

" ' yisiinpe~ ie e~ 2lv cos 93 y ' ^ ' 

Then the parametrization in Eq. (4) arises from an identification of E, P and D with g through the relation 

9' X dg + ig- x Tg - £T = ( jj. ) , ^Tg = -t (p. _0 (16) 

where we have used the constraint equation (12). Also, the zero curvature equation (10) with the identification in 
Eq.(16) becomes 



E* -ipX-iU lU \\P* -D 



= (17) 



whose components agree precisely with the SIT equation in the sharp line limit. The constraint equations (11) and 
(12), combined with Eq.(14), also reduce to Eqs.(7) and (8). Thus we have shown that the SIT equation consistently 
arises from the action (9) with the gauge fixing in Eq.(14). Moreover, the zero curvature equation (17) demonstrates 
the integrability of the SIT equation. The potential term in Eq.(9) changes into the population inversion D, 

/3 f j3 

-cos2v = J-D, (18) 

which for /3 = 1 possesses degenerate vacuua at 

cp = ip n = (n + tj) 71 "' n ^ ^ an< ^ ^ = ^0 for #0 constant . (19) 

The soliton solutions interpolating different vacuua can be obtained either by applying the dressing method or by 
using the Backhand transformation ||. In particular, the 1-soliton solution is given by 

COS(fi = — sec.Vi(2hz — 2bCz) 

rc = (a - + (a - £)Cz 

9 = - tan" 1 [^— ^-coth(2&z - 2bCz)] - 2£z + 2Dz (20) 



where a, b are arbitrary constants and 



In terms of E, 

E = -2ib sech(2bz - 2 bCz)e- 2z(az - Dg+{a -^ Cs) . (22) 

If a — £ = 0, this solution interpolates between two different vacuua ip n and <p n +i, i-e. it becomes a topological 

1- soliton (An = 1). On resonance where a = £ = 0, Eq.(20) reduces to the 1-soliton of the sine-Gordon equation, or a 

2- 7T pulse of SIT. If a — £ ^ 0, the solution in Eq.(20) reaches to the same vacuum asymptotically as x — > ±00 so that 
the topological number is zero (An = 0). Nevertheless, except for the topological number, this solution possesses all 
the properties of a soliton so that we call it a nontopological 1-soliton. It represents a localized pulse with a steadily 
varying phase. The time area of the pulse, which is defined by 

AS = 2 I \E\dt, (23) 



3 



is still 2tt. It is important to mention that this 2n area is a mere coincidence and should not be confused with the 
2-7T area of the topological one. Because of the interference between phases of each nontopological solitons, multi- 
nontopological solitons in general do not possess the area which is an integer multiple of 2ir and the area theorem of 
McCall and Hahn in the case of inhomogenous broadening does not hold. The stability of nontopological solitons, 
unlike the topological case whose stability is due to the topological protection, arises from the U(l)-charge conservation 
law. Recall that the action (9), consequently Eqs.(5)-(8), are invariant under the axial vector transform g — > hgh or, 
equivalently, 

X] — > r\ + e for e constant. (24) 

The corresponding Noether currents and the charge are given by 

cos 2 ^ - cos 2 ip = 

J = -7-5 — or/ , J = — orj 
sin ip sin ip 

/oo 
dx(J + J) (25) 
-OO 

where J, J satisfy the conservation law, dJ + dJ — 0. In particular, the charge of the 1-soliton in Eq.(20) is 

Ql-sol = -c(sign[6 ■ (a - 0]| - tan" 1 ^). (26) 

The stability of nontopological solitons can be shown either by using conservation laws in terms of charge and energy 
as given in || , or by studying the behavior against small fluctuations (jl0| . 

The action (9) also possesses two different types of discrete symmetries. This can be seen most easily in a different 
gauge where A = A = which is connected to the gauge in Eq.(14) by an appropriate vector gauge transformation. 
The first case is the chiral symmetry under the interchange, 

z <-> z and g «-> g~ x ( or r\ <-> —r\ , <p «-> —tp) (27) 

which is a characteristic of the WZW action. Unlike the case of a sigma model without the Wess-Zumino term, parity 
alone (z <-> z) is not a symmetry. In the SIT context, this is due to the slowly varying envelop approximation which 
breaks the parity invariance of the Maxwell-Bloch equation. The chiral symmetry generates a new solution from a 
known one. For example, the chiral transform of the 1-soliton in Eq.(20) in the resonant case (£ = 0) is again a 
1-soliton but with the replacement of constants a, b by 

(28) 



a 2 + b 2 ' a 2 + b 2 ' 

which changes the shape of the pulse as well as the velocity by v — > c — v. The currents and the charge also change 
into 

J — > — J , J — > — J , Q^-Q. (29) 
The second case is the dual symmetry of the Krammers-Wannier type || under 

<-> -0 , 9 i(J i9 ■ (30) 

Changing the sign of (3 makes the potential upside down so that vacuua become maxima of the potential and vice 
versa. This allows us to find a localized solution which approaches to a maximum of the potential asymptotically, i.e. 
it reaches to the completely population inverted state. In analogy with the nonlinear Schrodinger case, we call it a 
"dark soliton". In practice, the dark soliton for positive (3 can be obtained by replacing (3 — > —(3 , z — > — z in the 
(bright) soliton solution of the negative (3 case. For example, the dark 1-soliton can be written by 

cos^e 2 " 7 = b =tanh(2&z + 2bCz) - i , = 

V(a~0 2 + b 2 - 2 + b 2 

6 = -2{a-Ci{z-Cz)-2& (31) 

Finally, we show that inhomogeneous broadening can be incorporated into our formulation naturally with minor 
modifications. We maintain the constraint equation (12) only and modify the zero curvature eqution by 
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d + g- x dg + ig- x Tg - £T + AT , 5 



A-£ 



= 



(32) 



where the constant A is a spectral parameter which becomes A - 
an average over the spectrum /(£) as given by 



-£ in the sharp line limit. The angular bracket signifies 



< 



>= 



(33) 



We make the same identification as in Eq.(16) and require that g~ x dg + £,g~ x Tg — £T is ^-independent since E is 
a ^-independent macroscopic quantity. This results in the SIT equation with inhomogenous broadening as given in 
Eq.(l). Once again, by using the dressing method, the 1-soliton can be obtained which is the same as in Eq.(20) but 
with the replacement 



C = 



1 



(a - £) 2 + b 2 



D = (a-Z) 



1 



(a - 2 + & 2 



(a - O 2 + b 2 



(34) 



Eq.(14) shows that each frequency £ corresponds to a specific gauge choice in our formulation therefore inhomogenous 
broadening is equivalent to averaging over different gauge fixings. This implies that the inhomogenously broadened 
case can not be treated by a single field theory. Nevertheless it is remarkable that the group theoretical parametrization 
of E, P and D is still valid. Another important feature of inhomogenous broadening is that it introduces an anomaly 
term M in the U(l)-current conservation such that dJ + dJ — M and 



M = 2 cot ip[ cos(# - 2rj) < sm(6 - 2rj) sin 2ip > — sin(0 - 2rj) < cos(# - 2rj) sin 2tp > 
-(cot 2 tpdr) + -B6)d(p ]. 



(35) 



This anomaly vanishes in the sharp line limit due to the constraint Eq.(8). It also vanishes in the case of 1-soliton and 
the charge remains conserved. This may be compared with the conserved area of topological solitons in the presence 
of inhomogenous broadening. It is an open question whether there exists a similar theorem to the area theorem 
concerning about the stability of pulses with phase variation in terms of charge and anomaly. 

In conlusion, we have shown that the SIT equation in the sharp line limit can be identified with the complex 
sine-Gordon equation. This allowed us to present various mathematical properties of solutions as well as the equation 
of motion itself. Also, we have shown that the effective field theoretic description can be extended to inhomogenously 
broadened case by averaging over different gauge fixings. One of the advantage of the present reformulation of the SIT 
equation is that it clarifies the analytic nature of SIT in terms of group theroetical treatment and allows a systematic 
generalization to other multilevel cases than SIT. Also, the concept of U(l)-charge, which is a new feature of SIT 
coming from our for reformulation, can be used in explaining earlier numerical result of slow reshaping of pulses in 
the off-resonant case Q . Details will appear in a longer version of this Letter . 
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